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Abstract 

Four new exactly solvable, real and shape-invariant potentials as- 
sociated with a position-dependent effective mass are generated within 
the concept of shape-invariant potentials using a specific ansatz for su- 
perpotential. The accompanying energy spectra of the bound-state and 
the ground-state wavefunction are obtained algebraically as a function 
of free parameters and the results are compared with those of others 
works in the litterature. 
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1 Introduction 

Physical systems with position-dependent effective mass have received, in 
recent years, a significant attention due to their relevance in describing 
the physical properties of various microstructures such as compositionally 
graded crystals [1], quantum dots [2], semi conductor heterostructures [3], 
quantum liquids [4] and ^ff— clusters [5], etc. Recently, a wide number of 
exact solutions on these topics has increased [6-11]. In theoretical physics. 
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various methods and approaches are used inchiding Factorization and Oper- 
ator methods [12,13], Point canonical transformation methods [14], Super- 
symmetric quantum mechanical approach [15], Group theory approach [16] 
and Path integral formalism [17]. 

The effective mass Schrodinger equation is studied by many authors and 
its exact solutions (eigenvalues and eigenfunctions) are obtained. Quesne 
and Tkachuk have established a certain connection between the Schrodinger 
equation with position-dependent effective mass and the deforming func- 
tion appearing in the generalized canonical commutation relations. As a 
consequence, the potential in the deformed Schrodinger equation can be 
considered as the effective potentials in position-dependent effective mass 
[18,19]. 

The main purpose of the present paper is to obtain the energy spectra of 
the bound-states and the ground-state wavefunctions for four "new" solvable 
and real potentials (Three-dimensional Harmonic Oscillator, Morse, Poschl- 
Teller I and II, and three-dimensional Coulomb) deduced by applying the 
procedure of Ref. [18] in order to solve a position-dependent effective mass 
as a deformed shape-invariance condition introduced by Gendenshtein [20] 
and inspired by the supersymmetric quantum mechanical technique [21]. 
The one-dimensional potentials thus obtained are Shape-Invariant under 
parameter translations Aj+i = Aj it a where Aj and a denote two sets of 
three parameters, i.e. Aj = (Ai,(Jj,pj) and a = {a,b,c) with i = 0,1,2,... 
The new shape-invariance yields a new energy spectra of the bound-states 
with a non-equidistant spectrum which present some features with those 
obtained in Ref. [25]. 

The plan of the present paper is as follow. Factorization, effective poten- 
tial and Shape- Invariant Potentials are briefly reviewed in section 2. Section 
3 will deal with an assumption carried with the superpotential in order to 
discuss the shape invariance condition leading to the class of potentials and 
its corresponding "new" energy eigenvalues as well as their ground-state 
wavefunction. Finally, the last section contain the conclusion. 

2 Factorization, Effective potential and Shape In- 
variant Potentials. 

There are several ways to define the kinetic energy operator when the mass 
is a function of position. Since the mass and momentum operators are 
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no longer commute, the generalization of the standard Hamiltonian is not 
trivial. Therefore, the choice of the correct ordering of the operators of 
the kinetic energy to be Hermitian is indispensable. Defining a general 
Hermitian effective mass Hamiltonian proposed by von Roos [3] 



HvR = - (jn" {x)pm^ {x)prn^ {x) + rrC {x)pm^ {x)prrf (x)^ + V 



d 



(2.1) 



where /i = 1, a + /3 + 7 = — 1, and p = — i— — . The limits of the choice of the 

dx 

parameters a, (3 and 7 depend on the physical system. Using the restricted 
Hamiltonian from the a = 7 = constraint used by BenDaniel and Duke 
[22], we can write (2.1) as 



nvR = -d:,u^ (x)a^ + y(x), 

with C/^ ix) = ^ ^, , . Here we have used abbreviation = 

^ ' 1m(x) dx 

The identity found upon the commutation relation 



(2.2) 



brings the kinetic term in (2.2) to 



d.yWix) - yWix)d^ 
ld£_ix) 



(2.3) 



dxU^ix)d^ = d^^/U{x)U {x) y/U{x)d^ 
/— - 1 dccU (x) 



U{x) 



1 d^U (x) 



where the prime " ' " refers to the derivetive of U (x) with respect to x. 
The Hamiltonian (2.1) becomes [18,19] 

nvR = -(Vu~(^d^Vu(x)y + Veff{x), (2.5) 
where the effective potential Ve/f (x) is defined following (2.1) and (2.4) as 

V,ff{x) = V{x)+Vu{x), (2.6) 
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with 

V.W = ^"W^W+^. (2.7) 

In the formalism of supersymmetric quantum mechanics, there are two 
operators Q and Q^, called supercharges, that satisfy the anti commutation 
relation = Ti-sS' where T~Lss is the supersymmetric Hamiltonian 

[21]. The standard realization of the operators Q and Q'^ is Q = Aa- and 
Qt = where A (^^) and (T_ (cr+) are the bosonic operators and Pauli 

matrices, respectively. As a consequence of this, the appropriate operators 
to study the Hamiltonian (2.5) are 

Aeff = ^/ir(xjd^^irix) + Weff{x), (2.8.a) 
4// = -Vui^d^Vui^ + Weffix). (2.8.b) 

where Wg// (x) is called the effective superpotential. 

With this realization, the supersymmetric Hamiltonian (2.5) of the quan- 
tum system with position-dependent effective mass takes the form 

'Hl,eff = ^Iff^eff 

= -(ycy>)a,V[y>))' + yi,e//(a;), (2.9.a) 

where 

Vi,eff {x) = W^ff (x) - U {x) W^ff {x) , (2.9.b) 

and 

= -(^,/U{x)d^^/U{x)y + V2,eff{x), (2.10.a) 

with 

V2,eff {X) = W^ff {X) + U (X) W^ff (x) . (2.10.b) 

The Hamiltonian 7^2, e// is called the supersymmetric partner of Ti-i^eff ■ 
It can be easily shown that both TCi^eff Siud W2,e// has the same spectrum 
except for the ground-state, which belongs to Wi^e// [21]- It is obvious that 
(2.9.b) and (2.10.b) are related by 

V2,eff {X) = Fl,e// {X) + 2U {x) W^ff {x) . (2.11) 
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Substituting Vu {x) given by (2.7) as defined in (2.6) into Vi^eff (x) with 
i = 1,2 given by (2.11), we end up with 

V2ix) = Vi{x) + 2U{x)W^ff{x). (2.12) 

However, following the paper of Samani and Loran [23], the potential 
V2 (x) reads as 

V2 {x) = Vi (x) + 2U {x) W' {x) - U (x) U" (x) , (2.13) 

and comparing (2.12) to (2.13), we obtain after integration the relationship 
connecting both the superpotential and the effective superpotential 

W{x) = W,ff{x) + ^. (2.14) 

Inserting (2.14) into (2.9.b) and (2.10.b), we get [23] 

Vi{x) = W'^{x)-[U{x)W{x)]' , (2.15.a) 
V2 (x) = (x) - [U (x) W (x)]' + 2U (x) W' (x) 

-U{x)U"{x). (2.15.b) 

Despite their similar bound-state energy spectra, supersymmetric part- 
ner potentials constructed from (2.9-lO.b) and (2.15.a-b) usually have dif- 
ferent structures. However, the above potentials are called shape-invariant 
if V2 (x) has the same functional dependence on the coordinate as Vi (x) 
and differ only in some parameters [21]. The shape- invariant potentials are 
defined by the relationship 

V2 (x, ao) - Vi (x, ai) = 2U (x) W' (x, gq) - U (x) U" (x) 

= 7^(ao), (2.16) 

where ai = y (aq) is a function of parameters ao and TZ (ag) is independent 
of variable x. As a consequence, it can be shown that the discrete energy 
spectrum of Vi (x) can be written as [21] 

n— 1 

£;, = ^7^(afe), (2.17) 

fe=0 
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where the parameter is generated by the consecutive apphcation of the 
function y (x) , i.e. 

ofe = y^*"^ (oo) = yoyo---oy (00) • (2.18) 

^ V ' 

k-times 

Having determined the bound-state energy spectrum, the ground-state 
wavefunction of the corresponding Hamiltonian is obtainable by solving the 
first-order differential equation 

Aff% i^) = 0, (2.19) 

leading to the wavefunction 



Wn [xj = — , exp 



^ W(z) 
dz- 



U{z) 

where A/q is some normalization coefficient. 



(2.20) 



3 Exactly solvable potentials with three parame- 
ters 

The fairly general factorizable form of superpotential W {x) is given by Ref 
[21] 

W {x, ao) = 1] (fci + Ci) (x) + + h (x) , (3.1) 

1=1 

where ao = {ki, k2, . . .), ai = {ki + a, k2 + (3, . . .) and q, a, fJ, . . . being con- 
stants. It may be noted that one can obtain the condition to be fulfilled by 
the functions gi (x) , hi (x) and fi (x) once that (3.1) is used in (2.16). How- 
ever, Samani and Loran have carried out a simple form to W considering it 
in the ansatz with one-parameter [23] 

W{x,a) = ag{x) + ^+f{x), (3.2) 

where the functions g (x) , h (x) and / (x) are independent of the parameter 
a. 
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We shall now point out the specific ansatz carried out by the superpo- 
tential that goes into the determination of the effective potentials and its 
accompanying energy eigenvalues and ground-state wavefunctions. To this 
end, the "new" superpotential is characterized by one function and three 
parameters instead of that given by (3.2). To be more precise, we substitute 
all functions appearing in (3.2) by new parameters and the single parameter 
by a function, i.e. g (x) — > A, / (x) a, h (x) — > p and a ^ (p{x) , without 
making a point of seeking the transformation which connects them. 

Therefore the superpotential (3.2) becomes, taking into account (2.14) 

W (x. A) = A<^ (x) + ^ + a + (3.3) 

and the effective shape-invariance condition reads as 

V2,eff (x, Ao) = Vi,eff {x, Xi) + 7^ (oq) , (3.4) 

with Xi = {Xi,ai,Pi), i = 0,1,2,... Henceforth, we will suppose that A 
coincides with Aq. 

Inserting (3.3) into (3.4) we obtain 

(A? - Xl) 0^ (x) + 2 (Aiai - Aoao) 4> (x) + + 

(p (x) <p{x) 



U (x) (pi + TTpo) - U (x) (Ai + Ao) (x) 



<P'{x) 

4>\x) 

+ al-al + 2 (piAi - PqAq) + ^ (Aq) = 0. (3.5) 



with TT = ±1. The positive case (tt = 1) is what one ends up obtaining in (3.5) 
leading to the well-known results listed and tabulated in Refs. [18,19,21], 
while the negative case (tt = — 1) was added thinking that it leads to more 
interesting one. In the subsequent developments, we will be interested in 

the last case. 

Then, one way for (3.5) to be consistent and solvable is to separate the 
constant terms from the functions ; i.e. 



(A? - Xl) <t>^ (x) + 2 {X,a, - Aoao) 4> {x) + + + ^ 

^ (x) 0(X) 



U (x) (Ai + Ac) ct>' (x) - U (x) (pi - po) (3.6) 

4> {x) 
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and 

al-al + 2 (piAi - PqAq) + ^ (Aq) = q, (3.7) 

with q ^ 0. The resolution of (3.6) in terms of (/) (x) amounts to comparing 
these two members. Consequently, it is of primary importance to write a 
term U (x) 4>' (x) in such a way that the member of left-hand side is identified 
to that of right-hand side. To this end, three particular cases (solutions) 
arise, and henceforth, we will name them constant, linear and quadratic 
solutions and are identified with differential equations, respectively 

U (x) (p' {x) = a, (3. 8. a) 

Uix)(p'{x) = a(t>{x) + b, (a<0), (3.8.b) 
U{x)(l)'{x) = (x) + 60 (x) + c. (3.8.c) 

where a, b and c are constants. It is obvious that the solutions in (p (x) can 
be explicitly carried out by Euler's type integration [24]. 

3.1 Constant solution : Three-dimensional hcirmonic oscilla- 
tor 

3.1.1 Superpotential and Effective potential 

After a simple integration, the function cf) (x) is given by 

(pix) = an{x) + b, (3.9) 
where (a, b) G and jj, (x) is a function defined as a dimensionless integral 

z 

^^i-) = l^y (3.10) 

and which will appear frequently in subsequent subsections. 

Using both (3.3) and (2. 15. a) we obtain, respectively, the superpotential 
and the corresponding effective potential 

W (x. A) = A [ail (x) + b) + ^^^^^y—^ + <y + (3.11) 

V,ffix,X) = (a/x (x) +bf+ / + "I + 2Aa (a/x (x) + b) 

[afj. [x) + b) 

+ ^ + + 2Ap - aX. (3.12) 

a/i (x) + 
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It is obvious that if the a = b = constraint holds, then the effective 
potential (3.12) is related to shape-invariant three-dimensional harmonic 
oscillator potential 

yJH-O) {x, A) = A^aV' {x) + + 2Xp - aX. (3.13) 



a^//^ (x) 



3.1.2 Energy eigenvalues 



The energy eigenvalues can be calculated algebraically from (3.7). Indeed, 
inserting (3. 8. a) into (3.6) we obtain 

'^i - '^o = -^1^1 - -^oo-Q = Pi<^i - Po^o = 0, (3.14.a) 

Pi-Po = -a{pi-po), (3.14.b) 

q = a(Ai + Ao)7^0. (3.14.C) 

Solving (3.14) gets Pi = — (po + ) Ai = Aq and ai = ctq, which satisfies 
the recursion relations 



Pk 



{-'^f Po - |(l - (-1)'') > = Ao, = o-Q, (3.15) 



where k = 0, 1, 2, . . . Thus the energy eigenvalues are given, taking into 
consideration (2.17) 



n— 1 

£n = y^7^(Afc,crfc,Pfc) 
k=0 

n-1 

= ^4Afc(pfc + a) 

k=0 

n—1 n—1 

= 4a ^ Ao + 4^ Ao [(-1)^ - | (l - (-1)')" 

fc=0 fe=0 

= 2aAon + Ao(a + 2po)(l- (-1)"). (3.16) 

This spectrum presents some features with that obtained in formula 
(10) in Ref. [25]. Hence the negative case is allowed to generate new shape- 
invariant potentials with richer new energy spectra of the bound-states with 
a non-equidistant spectrum. 
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3.1.3 Ground-state wavefunction 



Since a particle is constrained to move in three-dimensional Harmonic Os- 
cillator, we set cTjfc = 6 = (A; = 0, 1, 2, . . .) , Aq = — and Pq = al {I <0) . 

Then, keeping in mind (3.10), the ground-state wavefunction is calculated 
straightforwardly using (2.19), we finally get 



Mo 



exp 



z z z 

\ f ^ /^(^) Po f dz 1 f 



Mo _i, . r w 2, X 
jj-^fi (x)exp[--// (x) 



dU{z) 
iHz) 

(3.17) 



3.2 Linear solution : Morse potential 
3.2.1 Superpotential and Effective potential 

In this case, the function (f) (x) is given by 

(p (x) = — (b — exp [—ap (x)]) . 



(3.18) 



with a ^ 0. The corresponding superpotential and the effective potential 
read as 



W (x, A) 



A 



exp [-ap (x)]) + 



ap 



b — exp [—ap (x) 



+ (3.19) 



Veff{x,\) = -^{b-exp[-ap{x)\) +- 

o [b — exp [—ap {x)\) 

ap(2a — a) . f 2a \ . , 

+ , ^ r f^ + A — + 1 6-exp -a/i X 

b — exp [—ap. (x)J \ a J 



+a'^ + 2Xp - bX. 



(3.20) 



From (3.20), if the p = b = constraint holds then the effective potential 
is related to shape-invariant Morse potential 

K/T'"^ (x. A) = ^ e-2«MW -A + 1^ e-'^'^(-) +a''. (3.21) 
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3.2.2 Energy eigenvalues 



X2 \2 



Inserting (3.8.b) into (3.6), we get the system of parametric equations ac- 
cording to 

0, (3.22.a) 

(3.22.b) 
(3.22.C) 
(3.22.d) 
(3.22.e) 



2 (AiCTi — Aqito) 

2 2 

Pi- Po 
2 - Po^^o) 



a(Ai + Ao), 
-biPi-po), 
-a{pi-po), 
q = 6(Ai + Ao)/0. 

and its solution lead to the following recursion relations 

Pk = {-'^f Po-\{}- i.-'^f) ^ Afe = Ao, ak = (TQ + ak. 
The energy eigenvalues can be calculated easily 

n-l 

= y^7^(Afc,crfc,Pfc) 
fc=0 
n-l 



(3.23) 



= ^ (46Afe - 2a(7jt + 46pfe - a^) 



fe=0 

= n(26Ao-2a(7o-a2n) +Ao(2po + 6)(l- (-1)") (3.24) 

Again, we obtain a richer spectra for the Morse potential as was done in 
(3.16) for the three-dimensional harmonic oscillator. The energy eigenvalues 
corresponding to the effective potential (3.21) can be obtained by setting the 
Po = ^ = constraint 

4Morse)^_^2(^^ + ^^%^2_ (325) 
3.2.3 Ground-state wavefunction 

The ground-state wavefunction ■i/^g {x) associated to the energy eigenvalues 
(3.21) is given after integration 



-00 [x) 



A/"o 
U{x) 



exp 



afi{z) 



a . U(z 



dz 



dU{z) 
lj{z) 



(3.26) 
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3.3 Quadratic solution : Poschl- Teller (I, II) and three di- 
mensional Coulomb potentials 

Integrating (3.8.c) in terms of the function (f){x), it is straightforward to 
obtain (sec 2.172 of Ref. [24]) 



I 



2 b + 2a6 (x) , „ x 

^ - : arctan ' , A > (3.27.a) 



«r + &C + c VA VA 



-2 , b + 2ad)(x) ^ ^ ^„ , x 
arctanh \ A < (3.27.b) 



-2 



- , ^ , / V , A = (3.27.C) 

b + 2a^{x) ' ^ ^ 

where A = — 6^ + 4ac. It is obvious that the integral appearing up is equal 
to a dimensionless mass integral, i.e. 

Wc will see in subsequent developments that the first two cases are re- 
lated to trigonometric and hyperbolic Poschl- Teller potentials, respectively, 
while the last case is associated to the three-dimensional Coulomb potential. 

3.3.1 Poschl-Teller I and II 

Superpotential and Effective potential The solution in terms of cf) (x) 



VA VA b 
is given through (3.27.a) by (f>{x) = — — tan— — /u(x) — — , the superpo- 

tential and the effective potential read as 

ly (X, A) = A f ^ tan (x) - ^ +^ ^ 

V- 2 2aJ #tan^Ma^)-^ ^ 

(3.28) 

A ^/A , , 6 y p2 



K// (a;, A) = A — tan —\x {x) - — \ + 



2 2a) (^tan#Mx)-A) 

/VA VA , , 6 \ 2po 
-\-2Xa \ — — tan— — /x(x) — 



2 



,2a 2 2a j ^tan^/.(x)-A 

AA pA 

4a _j_ 4a 



cos2^^(x) ' cos2^m(2:) (#tan^//(x)- 1^)' 
-Fcr^ 2Ap. (3.29) 



12 



The a = b = constraint leads to the shape-invariant trigonometric 
Poschl- Teller potential. Indeed, taking into account A = 4ac > 0, the 
effective potential is reduced to 

.(Trig.) , . _ ^ f, . ' 



cos^ yUici^ {x) sm"^ ^/ac|J, (x) \ \ a y cj 

The hyperbolic Poschl- Teller superpotential and effective potential are 
obtainable once the substitution -v/A i\/—A is made. As a consequence 
of this, the superpotential and the effective potential become respectively 

W{x,X) = -Af^^tanh^^/x(x) + ;^^ ^ 



2a 2 2a J ^ tanh (x) + A 

(3.31) 



Veff{x,X) = X' 



-^tanh^^(x) + - 



-2Xa 



2a 

AA 

4a 



tanh ■ 



(^tanh^/x(x) + Ay 
2po- 



2a 
pA 
4a 



tanh ■ 



_j_ ^±U 

cosh^ (x) cosh^ (x) tanh (x) + ^) ^ 



-Fcr^ 2Xp. 

As in the trigonometric case, the a = b = constraint leads to the shape- 
invariant hyperbolic Poschl-Teller potential. Taking into account A = 4ac < 
0, the hyperbolic effective potential is reduced to 



(3.32) 



V 



■^•'^ ' cosh^ ^J—acji (x) sinh^ ^—ac^ (x) 



Ac(^-l) 



- A 



- P 



(3.33) 
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Energy eigenvalues Substituting (3.8.c) in (3.5) we obtain the system 
of parametric equations 



^2 \2 



2 (AiCTi - AoCTo) 

2 2 

Pi -Po 
2 - poc^o) 



a(Ai + Ao) , 
6(Ai + Ao), 
-c(Pi - Po) , 
-b{pi - Po) , 

c(Ai + Ao)-a(pi-/)o)/0. 



(3.34.a) 
(3.34.b) 
(3.34.C) 
(3.34.d) 
(3.34.e) 



From (3.34), we deduce the following recursion relations 



Pk 



-1) Po 



1 



-1) 



Afc = Ao + ak, ak = ctq + bk. (3.35) 



then the energy eigenvalues are given by 
n-l 



^Tl{Xk,(Tk,Pk) 

k=0 
n-l 

J2 (4cAfc + 4apfc + AXkPk - "^bak + A) 



k=0 



= (apo + 2AoPo + cAo + y) (1 " ("1)") " ^ (2aPo + ac) (-1)" 
+n (2cAo + ac- 2bao) + (ac - }?) (3.36) 

Taking into account the constraint = 6 = as well as the sign of the 
discriminant the energy spectra of the bound-states for trigonometric 
and hyperbolic Poschl- Teller potentials are 



^•(Trig.) 



^(Hyp.) 



1 Ao ^ 
— V n ■ 

2 a 



aPo + 



+ 



apo 



(3.37) 



(3.38) 



The energy eigenvalues thus obtained present, as for the case of three- 
dimensional harmonic oscillator, some similar features with that obtained 
in formulas (22) and (29) in Ref. [25], respectively. 



^This means that sign a = signc for A > and sign a = — signc for A < 0. 
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Ground-state wavefunction The ground-state wavefunetion tJjQ (x) as- 
sociated to the trigonometric and hyperbohc Poschl- Teller potentials (3.30) 
and (3.33) are given respectively 

^ pQgA/a gj^-p/c _ (3 39) 

u (x) 

^(Hyp.) ^ cosh^/« (a;) sinh''/'^ v^A* (^) • (3-40) 

3.3.2 Three-dimensional Coulomb potential 

Superpotential and Effective potential Using (3.27.c), the function 
(pix) becomes 

and the corresponding superpotential and the effective potential read 

W (X, A) = -X ) - + <^ + (3.42) 

T/ ,2^ 2 + 5/i(.T) V , 4oVV(x) /2 + 6/i(.T) 



2afi{x) J {2 + bnix)y \ 2a/x (x) 

- 2TM;;) - + (2+5^(x))^ + + ^'-"'^ 

At first sight, the p = constraint leads to the shape- invariant three- 
dimensional Coulomb potential 

(3.44) 

Energy eigenvalues Since we are dealing with three-dimensional Coulomb 

1 A /A A , 2\ fbX ^ 

potential, we will assume that — l =t(t-|-l) and — a 

a \a J a \2a 

—Ze^, it is straightforward to get 

X = a{l + l), a = ^(/ + l) + ^^^, (3.45) 
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where I is the angular momentum quantum number, Z the atomic number 
and e the electronic charge. Here the parameters A and a coincide with 
Ao and (Tq, respectively. As a consequence to (3.45), the recursion relations 
(3.35) give 



Xk = a{k + l + l) , 



'^^ = 2(^ + ^) + 2(7Ti)+''- 



(3.46) 



for k = 0,1,2, .. . Since the potential (3.44) is obtained from the restriction 
A = pi^ = 0, this requires that the coefficients a, b and c are nonzero and pi^ = 



(-l)'po-^(l 



-11 



= imposes the condition that k should be even, 
i.e. k = 2p with p G N. As defined above in (3.36) and taking into account 
the last restriction, the shape- invariance condition function TZ {X^, crtj Pk) 
satisfies 



T^{\hCFkl) = 4:cXki-2baki 

= 4ac{k + l + l) -2b 

ZeH 



bil + l) 



+ 



Ze2 



2(Z + i; 



+ bk 



Vk + 



(3.47) 



The energy eigenvalues can now be obtained from equations (2.17) and 
(3.47). By inserting k = 2p m. (3.47), the eigenvalues can be rewritten as 



N-l 



k=0 



mt 



where int 
tion, int 



= - E 

p=0 

- b 

T+T 

N-r 



l + int 



,2 Ze^b 
N-l 



Zf? + 6 (Z + 1) int 



N-l 



3.48) 



is the greatest integer not larger than 



A^- 1 



-. By defini- 



X X — 1 

— for X even and for x odd. Here, both I and N are 

2 2 
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related by the relationship 



int 



"A^- 1" 



TV- 1 - s 



2n - 1 - (2/ + 1) 



n 
rir 



2 
1 



(3.49) 



with s = 21 + 1, n is the principal quantum number and a quantum 
number which denotes the number of radial nodes for the wavefunction. 
Here, we have assumed N = 2n such that the number of bound-state levels 



is equal to int 



1 

n + - 



= n for n 7^ 0. However, we can note, from (3.49), 



that the angular momentum quantum number I fulfills the condition I ^ 
1 3 

~2'~-^'"2''" 



2(n^ + Z + 1) 2(Z + 1) 



- , then the energy 



If wc make the replacement h 
eigenvalues (3.48) become 

- 4(„:r;l'imi)- + + + + ^ + " + • 

(3.50) 

For large values of I and n^, we can make the approximation Z ~ 1 and 
~ + 1 leading to impose that both quantum numbers take all integral 
values from to ?max and 

^rmax; respectively, i.e. \ — 0, 1, 2, ... , ^niax and 
Ur = 0, 1,2, . . . ,nrmax- Thus, the energy eigenvalues can now be obtained 
from (3.50) as 



c-(Cb.) 

''rirJ, 



Z'e 



2^4 



A{nr + l + lf {l + lf 



Ze^ - kF {rir, I) 



nl [ur + 21 + 2y 



2 (rir + l + l) 



with F {rir, = + + 1) + 1) and k 



(3.51) 



If! 
/ + 1 



. Prom general consider- 



ations, it is evident that the spectrum of negative eigenvalues of the energy 
will be discrete, while that of the positive eigenvalues will be continuous. It 
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follows from (3.51) that the function F {rir, I) is fulfilled by the condition 



F {rir, = + + 1) {2nr + 1) < . (3.52) 

K 

The energy spectra given in (3.51) have already been established by 
Quesne and Tkachuk [18,19]. 



Ground-state wavefunction Let us now complete this algebraic deter- 
mination of energy eigenvalues by a construction of the corresponding "ra- 
dial" ground-state wavefunction tpoi (x). Using the last considerations, the 

Tx. / X ^ + 1 Ze^ U' (x) 
superpotentiai can be rewritten as Wi [xj 



+ ■ 



■ + 



H{x) 2(Z + 1) ■ 2 
the "radial" ground-state wavefunction is given through (2.20) by 



and 







Uix) 



{x) exp 



2 (Z + 1) 



H{x) 



(3.53) 



4 Conclusion 



In the present paper, we have generated four "new" solvable, real and shape- 
invariant potentials simply by judicious applications of the supersymmetric 
quantum mechanics formalism and shape-invariant potentials. In all cases, 
we have derived the effective potentials as well as their accompanying en- 
ergy spectra of bound-states and ground-state wavefunctions. However, the 
new and the important contribution of this paper is to point out how the 
simple fact of going from the positive case (vr = 1), characterized by an 
equidistant spectra, to the negative one (tt = — 1), allows to generate a new 
shape-invariant three-dimensional harmonic oscillator, Morse and Poschl- 
TcUcr (I and II) potentials with non-equidistant spectra (3.16), (3.24) and 
(3.36), respectively, while the three-dimensional Coulomb potential has only 
a finite number of bound-states (3.51) in contrast with the standard coulomb 
problem. 

To conclude, the supersymmetric quantum mechanics can also be a useful 
machinery for the treatment of wide classes of potentials. It may be pos- 
sible to explore, in the context of position-dependent effective mass within 
the framework of non compact SO (2, 2) Lie algebra, all shape-invariant po- 
tentials listed and tabulated in Refs. [19,21], respectively. The works are in 
progress and will be deferred to later publication. 
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